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The spectra of soft gamma-ray repeaters (SGRs), with the excep- 
tion of the March 5, 1979 main burst, are characterized by high-energy 
cutoffs around 30 keV and low-energy turnovers that are much steeper 
. than a Wien spectrum. Baring (fil) found that the spectra of cas- 

' cades due to photon splitting in a very strong, homogeneous magnetic 

field can soften spectra and produce good fits to the soft spectra of 
SGRs. Magnetic field strengths somewhat above the QED critical field 
strength B CT , where B CI — 4.413 x 10 13 G, is required to produce cutoffs 
at 30-40 keV. We have improved upon this model by computing Monte 
Carlo photon splitting cascade spectra in a neutron star dipole mag- 
netic field, including effects of curved space-time in a Schwarzschild 
metric. We investigate spectra produced by photons emitted at dif- 
ferent locations and observer angles. We find that the general results 
of Baring hold for surface emission throughout most of the magneto- 

■ sphere, but that emission in equatorial regions can best reproduce the 
' constancy of SGR spectra observed from different bursts. 
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The association of supernova remnants with at least two of the three known 

■ SGRs (2-4) (SGR1806-20 and SGR0525-66, the Mar 5, 1979 source) is now 
Q-f well established. The third, SGR1900+14, is near a ROSAT source that is 
Q ■ possibly a supernova remnant (||). These associations strongly indicate that 

| this class of 7-ray bursts is linked to relatively young galactic neutron stars. 
^2 ' Furthermore, it has been suggested (^|) that SGRs are a special class of neutron 
stars, known as "magnetars", that have extremely strong magnetic fields. 
Those neutron stars that are born with periods of several ms can aquire high 
fields during a short period of vigorous convection that follows core collapse, 
j^j ■ The convection generates a rapid dynamo which can generate dipole fields as 
C3 ' high as 10 15-16 G. It has been pointed out that there are several attractive 
features of very high fields in accounting for observations of SGRs. Among 
these are (i) an explanation of the 8 sec periodicity of the Mar 5 event as dipole 
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spin down in the ~ 10 4 yr age of the N49 supernova remnant (ii) a reduction 
of the Compton scattering opacity below the cyclotron fundamental, allowing 
a photon flux ~ 10 4 above the Eddington limit (Q). 

Yet another advantage of very strong magnetic fields in models for SGRs (|l|) 
is the effectiveness of photon splitting in producing the soft spectra of SGRs. 
Photon splitting, 7 — > 77, attenuates 7-ray photons, degrading them to lower 
energies, where they split repeatedly until they escape the high-field region. 
Baring ([!]) showed that the emerging spectra of such photon splitting cascades 
could account for both the shape and the softness of observed SGR spectra if 
the field in the emission region B j£ 4 B cl . The splitting cascade spectra were 
computed by the solution of a kinetic equation for the photons, assuming 
a uniform field in a region of size R = 2 x 10 6 cm. This approach, while 
important in demonstrating the effectiveness of photon splitting in modelling 
SGR spectra, neglected the dipole field geometry and the strong gravitational 
field of a neutron star magnetosphere. We have improved upon the model of 
Baring ([!]) by including these effects in a Monte Carlo calculation of photon 
splitting cascades near strongly magnetized neutron stars. 



II. PHOTON SPLITTING CASCADE SPECTRA 



Photon splitting is forbidden in field-free regions but is allowed in neutron 
star magnetic fields. The splitting rate (^) for a photon of energy e in units 
of mc 2 , T sp (e) oc e 5 (B/B CI ) 6 sin 6 #kB, where #kB is the angle between the 
photon momentum and the magnetic field, can be large if B is near B cl -. 
The above rate is valid in the nondispersive limit, where £sin#kB ^ 2 and 
B <C B CI • It is dependent on the two polarization states of the photons in the 
bircfringent, magnetized vacuum: || or _L, where the photon's electric vector 
is parallel or perpendicular to k x B. The modes ||— >_L||, _L^_L± and -L— >|||| 
are the only ones that are non-zero through CP invariance. Since there is 
a difference in the rates of these modes which depends on field strength (|l|), 
polarized photons emerge from the emission region. In our calculations, we 
have included high field corrections to the above formula for splitting (^) that 
cause the attenuation coefficient to saturate somewhat above B ~ 4£> cr . 

Photon splitting cascade spectra will turn over roughly at the escape en- 
ergies of the photons, above which photons undergo at least one splitting 
generation, but below which the optical depth is always <C 1 and photons 
can escape the magnetosphere. The existence of such an escape energy is a 
consequence of the r~ 3 decay of the dipole field. Our previous analysis (Q) 
showed that escape energies for photons emitted at the neutron star surface 
and propagating through a dipole magnetic field are quite sensitive to the 
propagation angle 6^ of the photon to the dipole axis at the magnetic pole 
(6 = 0°), but nearly independent of this angle near the equator (9 — 90°). 
This effect is a function of the dipole field curvature: near the pole, the field 
lines are diverging rapidly, and different emission directions sample very dif- 
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FIG. 1. a) Ratio of the energies below which photons escape a dipole field without 
splitting in curved to flat space as a function of magnetic field, b) Photon splitting 
cascade spectra in curved and flat space for different surface field strengths. 9 and 
4> are the colatitude and azimuth angles of the photon emission point in the dipole 
field and # k is the photon emission angle to the dipole axis in the local inertial frame. 



ferent field orientations, but at the equator the field looks nearly the same in 
all directions. 
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In the present work we have included the general relativistic effects of curved 
spacetime in a Schwarzschild metric, following the treatment of Gonthier & 
Harding ( |i"o| ) who studied the effects of general relativity on photon attenu- 
ation via magnetic pair production. Those effects are the curved spacetime 
photon trajectories, the magnetic dipolc field in a Schwarzschild metric and 
the gravitational redshift of the photon energy as a function of distance from 
the neutron star. We have taken a neutron star mass, M = 1.4 Mq and radius, 
R = 10 6 cm in these calculations. The effects of curved spacetime decrease 
the escape energies by a factor of about 2 compared to flat spacetime (^), 
the largest contributions coming from an increase in the dipole field strength 
(by 1.4 at the pole) and the correction for the gravitational redshift, which 
increases the photon energy by roughly a factor of 1.2 in the local inertial 
frame at the neutron star surface. The qualitative behavior of escape energies 
as a function of 9^ in curved space is the same as in flat space: in the polar 
emission case, there is still a substantial variation in escape energy with emis- 
sion direction 9^ for both initial polarizations, but almost no variation in the 
case of equatorial emission. The ratio of curved to flat space escape energies 
as a function of magnetic field are shown in Figure la. At low field strengths, 
one is seeing the combined effects of the increase in the dipole field strength 
(which varies from pole to equator) and the increase in local frame photon 
energy; in high fields, only the increase in photon energy affects the escape 
energy, giving a factor of 1.2 decrease, due to the saturation of the splitting 
with increasing field strength. 

Figure lb shows Monte Carlo cascade spectra resulting from monoenergetic 
photons injected with e = 2 at the neutron star surface at the magnetic pole. 
The photons split many times before escape, each time dividing their energy 
into two photons with a distribution that peaks at half the parent photon 
energy, in the non-dispersive limit (8,9). The cascade spectra peak just below 
the escape energy for that field strength. Compared to those of Baring, these 
spectra show the same e 2 power law below the peak, but a more rapid decrease 
above the peak than the e -7 found by Baring, due to the fact that the Monte 
Carlo spatial injection occurs as a delta function at the surface, while the 
kinetic equation solution assumes an exponential injection. The peak in the 
flat space spectra in Fig. 1 are also a factor of around 2 higher than the 
homogeneous field case of Baring, due to the fall off in strength of the dipole 
field. The cascade saturates at B ~ 30£? cr , due to the saturation of the cross 
section in high fields, but the saturation energy as well as the peaks of all 
spectra are lower when curved space corrections are included. 

The effect of varying the emission location on the neutron star surface on 
the cascade spectrum is shown in Figure 2. Fig. 2a illustrates that the spectra 
are distinctly different for two emission angles at the magnetic pole, while they 
are almost identical for different observer angles at the equator. Thus, the 
behavior of the cascade spectrum is almost entirely determined by the escape 
energies at the initial photon injection point. 
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(a) 9 = 0, cp = 



(b) 9 = 90, cp = 90 
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FIG. 2. Photon splitting cascade spectra in curved space for photons emitted at 
the (a) pole and (b) equator of a dipole field and for different emission directions 9^ 
to the dipole axis. 



The generation of pure splitting cascades is contingent upon operation of 
at least two polarization modes of splitting so that polarization exchange is 
effected. Adler (||) demonstrated that in the dispersive magnetized vacuum, 
only one of the modes of splitting considered here, namely JL— ►||||, satisfies 
kinematic selection rules imposed by four-momentum conservation. This re- 
sult was derived assuming weak dispersion, i.e. a refractive index close to 
unity. While Shabad ([Tl]) has extensively looked at the regime of strong vac- 
uum dispersion, mostly near and above pair creation threshold, numerical 
computation of the refractive index and kinematic selection rules for super- 
critical fields well below pair threshold remains to be explored; this is a major 
goal of our future research. 

Should the selection rules extend to the strongly dispersive regime sampled 
by the SGR problem, cascade continuation could be effected by a polarization- 
switching process such as Compton scattering. The cross-section for this is 
suppressed below the cyclotron resonance (|l^ ) , and effective polarization state 
switching (||— >_L) will occur if the radiation is somewhat beamed along the field 
in the rest frame of the scattering electrons. Significant scattering opacities are 
quite plausible in the luminous environment of SGRs, and the complications 
of scattering in a dipole field geometry, which will tend to broaden the ouput 
cascade spectrum and smear its polarization spectrum, are deferred to future 
work. 

These calculations of photon splitting cascades in neutron star magneto- 
spheres lend support to the idea that the splitting mechanism in strong mag- 
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netic fields could be the cause for the softness and quasi-thermal shape of 
SGR spectra. The insensitivity of the cascade spectra to observer angle near 
the equator is very important for modeling SGR spectra. It implies that the 
observed spectra would not vary from burst to burst, even if the neutron 
star orientation changes (i.e. the star rotates), as long as the emission occurs 
at large magnetic colatitudes. This is consistent with the lack of observed 
spectral variation in bursts from SGR1806-20 fl3| ) and in the phase-resolved 
spectroscopy of the periodic soft tail of the Mar 5, 1979 event ([i"4|). 
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